Ballistic transport at room temperature in micrometer size multigraphene 
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The intrinsic values of the carriers mobility and density of the graphene layers inside graphite, the 
well known structure built on these layers in the Bernal stacking configuration, are not well known 
mainly because most of the research was done in rather bulk samples where lattice defects hide 
their intrinsic values. By measuring the electrical resistance through microfabricated constrictions 
in micrometer small graphite flakes of a few tens of nanometers thickness we studied the ballistic 
behavior of the carriers. We found that the carriers' mean free path is micrometer large with a 
mobility /i ~ 6 x 10 6 cm 2 /Vs and a carrier density n~7x 10 8 cm -2 per graphene layer at room 
temperature. These distinctive transport and ballistic properties have important implications for 
understanding the values obtained in single graphene and in graphite as well as for implementing 
this last in nanoelectronic devices. 

PACS numbers: 72.15.Lh, 72.80.Cw, 73.23.Ad, 81.05.Uw 
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The existence of massless-like quasiparticles (Dirac 
fermions) found in the quasi- two dimensional graphite 
structure^— and in single graphene layers^ enhanced 
enormously the interest in these materials, implying not 
only the start of a qualitative new physics^ but also the 
possibility of nanotechnological improvements 7,8 . The 
search for highest carrier mobility together with lowest 
carrier density in graphene remains a main research is- 
sue of this two-dimensional model system also because of 
their implications in future applications in nanoelectron- 
ics. It is nowadays known, however, that graphene prop- 
erties can be sensitively affected by the environment, par- 
ticularly by the substrate 9 -*^ and by disorder—. This fact 
is supported by experiments done on suspended graphene 
layers where mobility approaching ~ 5 x 10 5 cm 2 /Vs at 
carrier densities n ~ 10 9 cm -2 at T < 30 K has been 
reache d 12 i 13 . Nevertheless, the expectations on the use 
of single graphene layer (SGL) on substrates are high 7 
and in spite of the detrimental substrate influence, SGL 
exhibits ballistic transport in the 100 nm range at room 
temperature^ and in the ~ 1 /im range at very low 
temperature: 
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By measuring the resistance and its ballistic character- 
istics we studied the transport behavior of the carriers of 
the graphene layers inside the graphite structure. This 
study allows us to obtain the carriers' values for the mean 
free path, carrier density and mobility. Due to the natu- 
ral packing of the graphene layers inside graphite we ex- 
pect not only a natural shielding from the environmental 
influence but their should be actually more ideal than 
SGL standing alone or on a substrate. Although some 
experiments have been done in graphene bilayers, the 
actual intrinsic values of these transport parameters for 
the graphite structure remain still not well known simply 
because their are very sensitive to lattice defects^, grain 
boundarie o 17 ! 18 and internal interfaces^ 9 -, which in gen- 
eral trigger high carrier density with low mobility. Recent 
results obtained on bulk oriented pyrolytic graphite sam- 



ples indicate ballistic transport in the micrometer range 
at T < 10 K with a very low carrier density n < 10 9 cm~ 2 
and huge mobility \i > 10 6 cm 2 /Vsi£. 

The graphite flakes we report here were obtained by 
exfoliation of a highly oriented pyrolytic graphite sam- 
ple of ZYA grade. Using an ultrasonic technique we ob- 
tained several flakes that we selected by measuring the re- 
sistance, its temperature dependence and micro-Raman 
signals and taking into account the overall shape and di- 
mension of the sample, its surface and thickness. The Au 
contacts for longitudinal resistance measurements were 
prepared using electron beam lithography, see Fig. 1. 

The measurement of the carrier mean free path I is 
not straightforward and in general the value obtained 
depends on several others, not well known parameters 
within the selected transport model, usually based on 
a Drude-Boltzmann approach. However, in case of bal- 
listic conduction there is a transparent method to ob- 
tain £ without adjustable parameters, successfully used 
for macroscopic graphite samples as described inl&. The 
method is based on the measurement of the longitudinal 
resistance as a function of the geometry of a constriction 
located between the voltage electrodes (see Fig. 1). We 
prepared these constrictions with the focused ion beam 
of a dual beam microscope. We avoided the modification 
of the crystalline structure of the samples due to the ion 
beam spread by protecting them with a negative elec- 
tron beam resist (AR-N 7500) of ~ 300 nm thickness, 
a method successfully tested in graphite and described 

As shown iniS the resistance of a graphite sample of 
width O, thickness t, with a constriction of size W and 
length L connecting two half-parts of resistivity p is given 
by: 
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At the rhs of Eq.(l) the first term corresponds to the 
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FIG. 1. Scanning electron microscope images of the two sam- 
ples discussed here with the four electrodes and with one of 
the constrictions in between. The graphite flakes (white lines 
indicate their perimeter) cannot be seen directly because of 
the top 300 nm resin layer used to shield them from the Ga + 
focused ion beam. Two upper pictures correspond to sample 
A and bottom to sample B. 
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ballistic Knudsen-Sharvin resistance^; the second, loga- 
rithmic term to the ohmic, spreading resistance in two- 
dimensions with the smooth function j(k = W/i) ~ 
1 - 0.33/cosh(0.lK) = 0.67... 1 for k = 0...oc 21 . The 
last term is due to the ohmic resistance of the constric- 
tion tube itself. The constant a takes care of the influ- 
ence of the sample shape, the topology and the location 
of the electrodes in the sample. It is equal to 1 for the 
geometry used in sample A where the voltage electrodes 
are deposited through the whole sample width and their 
distance is ~ £1 ^> W. This can be verified from the 
measurement of the resistance as a function of W, as it 
is shown below. To keep the third term, the ohmic contri- 
bution of the constriction, small enough we prepared con- 
strictions with length L < 0.4 /im. The results presented 
here correspond to two samples, sample A and B with dif- 
ferent geometry and resistivity. Sample A (B) had a size 
(distance between voltage electrodes x width x thick- 
ness) of 12.7 x 16 x 0.015 /im 3 (2.7 x 9.2 x 0.040 /im 3 ) 
and a resistivity /?(300 K)= 89 /ificm (18 /iSlcm), see 
Fig. 1. 

The simplest and direct way to check the ballistic con- 
tribution and obtain the mean free path (without further 
adjustable parameter) is just measuring the resistance R 
as a function of the constriction width in combination 
with Eq. (1). The results for sample A at two tempera- 
tures are shown in Fig. 2. In that figure one recognizes 
that for W < 2 fj,m the ballistic contribution overwhelms 
the ohmic ones, see Eq.(l), indicating that the mean free 



FIG. 2. Measured resistance for sample A as a function of 
the constriction width W at (a) 60 K and (b) 250 K. The 
point with the largest W corresponds to the virgin sample 
without a constriction. The different lines correspond to the 
first ballistic term (dashed line) and the ohmic (2 (dotted 
line) and 3 r (dashed-dotted line) contributions in Eq. (1)) 
and the continuous line to the addition of the three. The 
only free parameter is the mean free path £, see Eq. (1). For 

(a) the continuous line is calculated with £ = 1.2 fj,m and for 

(b) I = 0.8 /im. 



path should be of this order. The theoretical lines in 
Fig. 2 were obtained using £ — 1.2 /im and 0.8 /jm at 
60 K and 250 K, respectively. We estimate an error of 
~ 35 % for these values determined by the errors in the 
measurement of the sample and constriction geometry. 
We can calculate now the Fermi wavelength per graphene 
layer using the relation 

X F = 2Tre 2 N s p£/ht (2) 

with N s the number of graphene layers in the sample. 
For sample A we obtain then = 0.5 (0.8) ± 0.25 ^m 
at 250 K (60 K). 

In case Xf is larger than the constriction size, the bal- 
listic contribution to the resistivity is better described 
by the inverse of a sum of an energy E- and transverse 
wave vectors qvj-dependent transmission probabilities T n , 
where n = 0, ±1, ±2, . . The maximum possible n is 
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determined by the constriction width W, decreasing the 
smaller W. In this case the increase in resistance is ex- 
pected to show an oscillatory behavior as a function of 
W or \p 23 i 24 as observed experimentally in Bismuth (Bi) 
nanowires^ as well as in GaAs device o 26 i 27 . 

Note that the distance between voltage electrodes in 
sample A is larger than the obtained t. The larger the 
sample the larger is the probability to have defective re- 
gions with larger carrier concentration and smaller mean 
free path within the voltage electrodes^. Therefore we 
repeated the experiment with sample B that shows lower 
resistivity and with a smaller voltage electrode distance, 
see Fig. 1. Figure 3 shows the measured resistance nor- 
malized by its value at a constriction W = 3 /xm; in this 
way we pay attention to the huge relative increase de- 
creasing W and we need neither the absolute value of p 
nor of a to compare the data with theory. 

In Fig. 3 one realizes that for sample B Eq. (1) does not 
describe the experimental data even assuming the largest 
possible mean free path equal to the voltage electrode 
distance of 2.7 /xm. The data can be reasonably well 
fitted dividing the ballistic term in Eq. (1) by the function 
tmnc(2W/ Xf)Xf/2W, which generates steps whenever 
the constriction width W ~ iX/2 with i an integer. From 
the fits we obtain the parameters Xf = 1-0 (1.5) /xm 
and I = 2.2 (2.7) ± 0.3 /xm at 300 (10) K, see Fig. 3. 
Using other values of £, for example I = 1.3 /xm, see 
Fig. 3, the function does not fit the data indicating indeed 
that the carriers behave ballistically between the voltage 
electrodes, leaving actually Xf the only free parameter. 

The ballistic analytical function we use resembles the 
theoretical results with similar steps obtained ir>2i where 
the conductance vs. W was calculated numerically for 
a SLG with an electrostatically potential landscape that 
resembles a constriction. An analytical average value or 
envelope of this stepped function is obtained replacing 
the truncation function by exp(— Xp/2W), see Fig. 3. 
This exponential function represents the impossibility 
of an electron to propagate in the constriction when 
W < Air/2, as also occurs for the propagation of light 
in a tube. 

The data appear to be better fitted by the stepped 
function than with the exponential one. The important 
result obtained for sample B is the huge increase of the 
resistance for W < 2 (im indicating clearly a larger i. The 
larger I in sample B to the one obtained in sample A is 
compatible with the measured resistivity difference. 

As a further proof that the huge increase of the re- 
sistance decreasing the constriction width in sample B 
is due to the ballistic contribution and not due to, e.g. 
a possible disorder produced by the ion beam on the 
graphite structure, we show in Fig. 4 the temperature 
dependence of the measured resistance without and with 
the different constriction widths. We observe that in 
spite of the huge resistance increase the temperature 
dependence remains similar. The small differences in 
the temperature dependence of the resistance, decreas- 
ing slightly more at T > 80 K at large values of W than 
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FIG. 3. Normalized resistance for sample B vs. constriction 
width W; note the semilogarithmic scale, (a) Data at 10 K. 
The line with steps is obtained dividing the ballistic term 
in Eq. (1) by (A F /2W)trunc(2M//A F ) with parameters I = 
2.7 /xm and Af ~ 1-5 /xm. The continuous line is obtained 
multiplying the ballistic term in Eq. (1) by the exponential 
function exp(AF/2W). The dotted line follows Eq. (1) with 
£ = 2.7 /xm. (b) The same as in (a) but at 300 K. The used 
parameters are Af = 1-0 /xm and £ = 2.2 /xm. The dotted 
stepped function is obtained using the same Af but with a 
smaller £ = 1.3 /xm. 



at the smallest W, e.g. R(2)/R(300) ~ 1.6 (sample with- 
out constrictions) to ~ 1.45 for W = 0.6 /xm (see Fig. 4) 
can be explained taking into account that the ballistic 
contribution (first term in the r.h.s of Eq.(l)) is weakly 
T dependent. At low values of W the ohmic contribu- 
tions start to be negligible compared with the ballistic 
one (see Fig. 2) and therefore the temperature depen- 
dence of R(T) slightly reduces decreasing W, see Fig. 4. 

As shown ir>i£ the temperature dependence of R(T, W) 
can be used now to obtain Xf{T) and the mobil- 
ity per graphene layer, this last given by /x(T) = 
(e/h)XF(T)£(T). Since the density of carriers per 
graphene layer can be calculated from n = 2ir/X F we 
show in Fig. 5 the mobility vs. carrier density for the two 
samples and compare them with data from literature for 
suspended SLG. From these results we clearly recognize 
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FIG. 4. Resistance of sample B vs. temperature at different 
constriction widths and without constriction ((+),f2 is the 
total sample width, i.e. data of the sample without constric- 
tion). 



the much larger mobility and smaller density of carrier 
for the thin graphite flakes, supporting the view that the 
graphene layers within graphite are of better quality than 
SLG and with a smaller carrier density. Obviously sam- 
ple B is not free from defects and therefore we expect 
that the obtained values might still be improved in ideal, 
defect free graphite structures. 

The rather weak temperature dependence of n(T) (or 
Xf) obtained for the here reported samples is a non- 
trivial result that deserves a comment. For the simple 
case of a linear dispersion relation we expect that the 
A F (T) oc 1/T and therefore n(T) oc T 2 or oc T for a 
quadratic dispersion relation. The data obtained indi- 
cate that the carrier density increases at most a factor 
of ten between 10 K and 300 K. We may speculate the 
reason for this behavior, either by comparing our results 
with those found in Bi and a possible band splitting or 
the existence of an intrinsic energy gap in graphite. 

Firstly, we note that both materials, graphite and Bi, 
have very small carrier densities with small anisotropic 
effective carrier mass. For a long time it was unclear 
why there was such a difference in the carriers' mean free 
path between the two structures. While the carriers in 
Bi show £(300 K) > 2 /xm ; 28 ' 29 the values reported for 
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FIG. 5. Mobility vs. electron density obtained for samples 
A and B. The points are obtained between 300 K and 10 K 
(60 K for sample A). The line corresponds to the data of a 
suspended SLG at 20 K from Ref. Q3. 



graphite in literature are more than one order of mag- 
nitude smaller at the same temperature. The results of 
this work finally put an end to this apparent difference 
showing that also in graphite the carriers' mean free path 
can be micrometer large with a low carrier density. 

The large Xf of the order of the sample dimensions 
may affect the band structure. For Xf ~ 1 /im, and for 
the typical sample dimensions we report here only the 
first five quantized wave vectors in reciprocal space are 
occupied (Ep ~ 4 meV) indicating that the usual con- 
tinuous approximation of the band structure is inappro- 
priate. As in Bi nanowires we expect that quantum con- 
finement affects the electronic band structure and a split 
into subbands can occur—. It is interesting to note that 
the temperature dependence of the resistivity of our sam- 
ples, see Fig. 4, is basically the same as for Bi nanowires 
of diameters <fi < 70 nrr>2£. Therefore, it may be possible 
that the carriers in the micrometer small graphite sam- 
ples fill partially a narrow band with energy gap of the 
order of ~ 10 meV, preventing a simple thermally acti- 
vated excitation of the carriers to the conduction band. 

Because macroscopic graphite samples have a much 
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larger density defects (n(T » 1 K) > 10 11 cm" per 
graphene layer) than our samples^, it appears plausi- 
ble that the temperature behavior of the resistance ob- 
served here cannot be seen in bulk HOPG. Experimental 
evidence for the influence of interfaces to the transport 
properties of HOPG was reported recently^ supporting 
this statement. It remains unclear what we expect to 
happen in ideal graphene. The reported data in liter- 
ature were taken mostly for samples with larger carrier 
density than in our samples. Therefore, we believe that 
one cannot clearly answer yet on the existence or not 
of an energy gap in ideal graphene (experiments at car- 
rier concentration < 10 8 cm~ 2 were not yet reported). 
At such low carrier densities, as it appears to manifest 
in graphite, electron correlations and possible localiza- 
tion effects should be considered. Electron interactions 
are large and for a small enough carrier density, the ex- 
pected screening will be very weak promoting therefore 
the existence of an energy gap. This is what it is observed 
in Monte Carlo simulations for the unscreened Coulomb 
interaction in graphene with Dirac flavor N s = 2, see 
Ref. M 

The obtained mobility in our thin micrometer small 
graphite samples is two orders of magnitude larger, with 
a carrier density three orders of magnitude smaller than 
those measured in bulk samples, see e.g. Ref. l32l indi- 
cating that macroscopic samples do not provide neces- 
sarily the intrinsic properties of graphite due in part to 
the large density of lattice defects^. The ballistic be- 
havior obtained here with micrometer large mean free 
path and low carrier density was further verified by mag- 
netoresistance and Hall effect measurements and by the 
observation of Aharonov-Bohm oscillations in the mag- 
netoresistance. 

The magnetoresistance for fields normal to the 
graphene planes systematically decreases decreasing the 
constriction width, starting already with the largest W 
as shown in Fig. 6 for sample A (similar behavior is ob- 
tained for sample B), as observed in bulk HOPG sam- 
ples but at lower temperatures^. This behavior indi- 
cates that I > 1 \xm up to room temperature. The field 
dependence of the Hall resistance as well as its temper- 
ature dependence in similar graphite flakes as presented 
in this work, deviate clearly from measurements in bulk 
graphite samples, see e.g. Ref. [H. Although the cal- 
culation of the carrier density through the Hall effect in 
semimetals due to the electron and hole contributions is 
non simple, the obtained results support not only the 
existence of the low carrier density but also the weak 
temperature dependence of n. Finally, for fields parallel 



to the graphene planes and input current the magnetore- 
sistance shows pronounced Aharonov-Bohm oscillations, 
similar to that found in Bi nanowiresi 4 - and topological 
insulators 3 ^, demonstrating the coherent propagation of 
electrons around a length of several micrometers up to 
temperatures T ~ 250 KM 
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FIG. 6. Magnetoresistance vs. magnetic field applied nor- 
mal to the main surface of sample A for different constriction 
widths and at constant temperature of 250 K. The clear re- 
duction of the magnetoresistance decreasing W supports the 
existence of a mean free path of the order of several microns. 
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